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ABSTRACT 

We present here the first convincing observational manifestation of a magnetar-like 
magnetic field in an accreting neutron star in binary system - the first pulsating ultra- 
luminous X-ray source X-2 in the galaxy M82. Using the Chandra X-ray observatory 
data we show that the source exhibit the bimodal distribution of the luminosity with 
two well-defined peaks separated by a factor of 40. This behaviour can be interpreted 
as the action of the “propeller regime” of accretion. The onset of the propeller in a 
1.37 s pulsar at luminosity of 10“^° erg s“^ implies the dipole component of the 
neutron star magnetic field of ^ 10^"^ G. 

Key words: accretion, accretion discs - magnetic fields - stars: individual: M82 X-2 
- stars: magnetars -- X-rays: binaries. 


1 INTRODUCTION 

The revolutionary discovery of the pulsating ultra-luminous 
X-ray source (ULX) X-2 in the galaxy M82 (known also as 
X42.3-1-59; Kaaret et al. 2006) made by the NuSTAR obser¬ 
vatory (Bachetti et al. 2014) brought more questions than 
answers on the physics of accretion onto magnetized neu¬ 
tron stars as well as on the nature of ULXs. The chief dis¬ 
tinction of this source is an extremely high luminosity for an 
accreting neutron star of about 10'*° erg s“*. Other param¬ 
eters of the system, such as the spin period of P = 1.37 s, 
the orbital period of 2.5 d and the companion star of 5.2 Mq 
(Bachetti et al. 2014), are quite typical for an accretion pow¬ 
ered X-ray pulsar in a high-mass X-ray binary. The variable 
over the NuSTAR observations spin-up was also detected 
with an averaged rate P ~ —2 x 10~*° s s“*. 

Exceeding the Eddington luminosity by almost two or¬ 
ders of magnitude is a challenge for current theories and 
imposes strong limitations on the physical conditions in the 
emitting region. Many different models are already proposed 
aiming at the understanding of the nature of this source. 
These models can be divided into two main groups depend¬ 
ing on the neutron star magnetic field. One of the possi¬ 
bilities to emit such a high flux is to assume a magnetar- 
like magnetic field of the order of 10*"* G to significantly 

* E-mail: sergey.tsygankov@utu.fi 


reduce the interaction cross-section between radiation and 
the infalling material (see, e.g., Harding & Lai 2006). How¬ 
ever, such a magnetic field is almost two orders of magni¬ 
tude higher than the typical value observed in X-ray pulsars 
(Walter et al. 2015) and cannot be probed through spectral 
analysis. Therefore, all conclusions made so far are done 
based on the timing properties of M82 X-2. Particularly, 
assuming torque equilibrium, Ek§i et al. (2015) got the B- 
field strength of (2 A 7) x 10*° G. Dall’Osso et al. (2015) 
favoured a lower magnetic field of B ~ 10*° G, solving the 
torque equation numerically. 

At the same time other authors declare a low magnetic 
field in M82 X-2 based on the same timing properties of the 
source. For instance, Kluzniak & Lasota (2015) argue that 
the observed torque is consistent with the accretion disc ex¬ 
tending down to the vicinity of the neutron star surface. This 
can be the case only if the dipole magnetic field of the star 
is low, B ^ 10° G. In the original paper by Bachetti et al. 
(2014), the magnetic field of the order of 10*° G was esti¬ 
mated under the assumptions of the spin equilibrium and 
the Eddington accretion rate. 

Thus, the lack of an obvious observational manifesta¬ 
tion of the strong magnetic field precludes from any final 
conclusion. Here, using the Chandra observatory data, we 
show that the ULX M82 X-2 regularly enters the “propeller 
regime” of accretion which is seen as dramatic variations of 
the emitted luminosity. These observations imply the dipole 
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Table 1. Chandra observations of ULX M82 X-2. 


Obs Id 

Instrument 

Date 

(MJD) 

Exposure 

(ks) 

Pile-up fraction 
(per cent) 

a 

Photon 

index 

Luminosity^ 
(10^® erg s-l) 

361 

ACIS-I 

51441.474 

33.25 

1 

- 

2.1 ±0.4 

2.6 ±0.5 

1302 

ACIS-I 

51441.883 

15.52 

1 

- 

1.3 ±0.7 

2.0 ±0.6 

1411_000 

HRC 

51479.184 

36.04 

- 

- 

0.4 ± 1.2 

114.3 ±11.3 

1411_002 

HRC 

51563.619 

17.61 

- 

- 

1.35*' 

1.7± 1.0 

2933 

ACIS-S 

52443.784 

18.03 

20 

0.24 ±0.21 

1.4 ±0.2 

60.9 ±5.4 

5644 

ACIS-S 

53599.038 

68.14 

5 

- 

1.4 ±0.1 

107.4 ±3.4 

6361 

ACIS-S 

53600.664 

17.45 

5 

- 

1.3 ±0.1 

105.1 ± 3.3 

8189 

HRC 

54109.345 

61.29 

- 

- 

l.SS*” 

3.2 ± 1.7 

8505 

HRC 

54112.093 

83.22 

- 

- 

l.SS*” 

3.4 ± 1.5 

10542 

ACIS-S 

55006.163 

118.61 

1 

- 

1.0 ±0.2 

3.8 ±0.4 

10543 

ACIS-S 

55013.936 

118.45 

1 

- 

0.7 ±0.4 

2.3 ±0.4 

11104 

ACIS-S 

55364.136 

9.92 

30 

0.39 ±0.09 

1.1 ±0.2 

119.5 ± 9.6 

13796 

ACIS-S 

56148.648 

19.81 

30 

0.23 ±0.11 

1.1 ±0.1 

173.8 ±18.4 

15616 

ACIS-S 

56347.964 

2.04 

5 

- 

1.35*’ 

6.4 ±2.3 

16580 

ACIS-S 

56691.841 

46.85 

30 

0.43 ± 0.23 

1.2 ±0.1 

187.4 ±24.4 


“ Luminosity in the 0.5-10 keV energy range corrected for interstellar absorption and assuming 
distance to the source D = 3.3 Mpc (Foley et al. 2014). 

^ Fixed at the averaged value obtained in the bright observations least affected by the pile-up effect 
(Obs ID 5644 and 6361). 


component of the neutron star magnetic field of ~ 10^"^ G, 
independently confirming the magnetar nature of this ULX. 


2 CHANDRA OBSERVATIONS AND RESULTS 

The Chandra X-ray observatory monitored the galaxy M82 
more or less evenly during the past 15 years resulting in 29 
publicly available observations. Some of them were pointed 
quite far from the ultra-luminous X-ray source M82 X-2. For 
our analysis, we selected only on-axis observations where the 
PSF shape allowed us to confidently separate the flux from 
M82 X-2 from the nearby sources. To exclude the selection 
bias we checked that none of the observations was triggered 
on a particular state of this source. The final sample of the 
utilized observations consisting of 15 pointings is listed in 
Table 1, where observation Id, name of the instrument, date 
and exposure of observation, as well as pile-up fraction are 
given. The level of pile-up was estimated using the CIAO 
tool PILEUP_MAP, which calculates an average number of 
counts per ACTS frame. These numbers were then used to 
estimate the pile-up fraction.^ 

The reduction of the data has been done following the 
standard pipeline in the Chandra Interactive Analysis of Ob¬ 
servations software package (CIAO, version 4.7). The source 
flux was estimated from the energy spectrum extracted us¬ 
ing CIAO SPECEXTRACT tool from a circular region of 
radius ~1". The background was extracted from an annulus 
region centered at the source position with the inner and 
outer radii l'.'8 and T'!Q, correspondingly. All contaminating 
point sources were excluded. 

Spectral fitting was performed with XSPEC v. 12.8.Ig 
(Arnaud 1996) using a simple power-law model with inter¬ 
stellar absorption (phabs model). Spectra were groupped to 
have a minimum of 1 count per bin and the Cash statistic 

^ http://cxc.harvard.edu/ciao/download/doc/pileup_abc.pdf 


was used. Some observations in the high-luminosity state 
suffer from the effect of pile-up with a pile-up fraction 
around 30 per cent at worse. This makes a reliable spec¬ 
tral fitting difficult, also taking into account the photon in¬ 
dex and absorption values degeneracy in the Chandra energy 
band. The same is valid for low-luminosity states where we 
do not have enough counts to constrain spectral parameters 
independently. To solve this problem, we have selected ob¬ 
servations in bright states which are much less affected by 
pile-up (Obs ID 5644 and 6361) and obtained the value of 
the hydrogen column density. Ah = 3.11 x 10^^ cm“^, and 
used this value to fit all remaining spectra allowing only the 
photon index and normalization to vary. Such an approach 
is supported by the recent work by Brightman et al. (2015) 
where a detailed spectral study of M82 X-2 was performed 
and no significant variability of Ah and F values between 
different Chandra observations was found. Additionally we 
checked the consistency of the spectral shape in different 
intensity states by the simultaneous fitting of ali spectra ob¬ 
tained in the low state. The resulting spectral parameters 
are in very good agreement with ones in the bright state. 
For fitting the piied up spectra from observations with the 
pile-up fraction higher than 5% we followed the Chandra 
ABC guide to pile-up,^ namelly we added the pileup model 
(see also Brightman et al. 2015). The best-fit a parameter, 
characterizing the grade morphing, as well as the photon 
index are shown in the Table 1. 

The flux value and its la uncertainty were calculated 
using CFLUX model from the XSPEC package with the fixed 
hydrogen column density value. The resulting luminosity in 
the energy range 0.5-10 keV corrected for the absorption is 
shown in the Table 1. The luminosities of M82 X-2 measured 
by us are in a reasonable agreement with the recent results 
by Brightman et al. (2015), especially if one takes into ac¬ 
count the density of point sources and the presence of the 
diffuse emission in the region around X-2 (see Fig. 1). 

In order to estimate the bolometric ffux, we assumed 
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Figure 1. Chandra images of M82 galaxy’s centre during observations performed on September 20, 1999 (MJD 51441.47) when M82 
X-2 was in a low-luminosity state (left) and August 17, 2005 (MJD 53599.04) when it was in a high-luminosity state (right). Circles 
indicate the positions of M82 X-1 and X-2 ultra-luminous X-ray sources. 


the spectral shape of M82 X-2 to be typical to that of 
X-ray pulsars (power law modified by a high-energy cut¬ 
off at ~ 15 keV with folding energy ~ 15 keV; see, e.g., 
Filippova et al. 2005), resulting in a bolometric correction 
factor of 2. This factor is consistent with the broadband 
spectrum of the pulsed emission from M82 X-2 as seen by 
NuSTAR (Brightman et al. 2015). 

The light curve of M82 X-2 as observed by Chan¬ 
dra is shown in Fig. 2(a). The histogram of the luminosi¬ 
ties shown in Fig. 2(b) clearly demonstrates a bimodality, 
with two well defined peaks at ~ 1.0 x 10"^° erg s“^ and 
~ 2.8 X 10^® ergs~^.^ We stress here, that in the majority 
of low-luminosity states the source is still presented in the 
Chandra images. This can be illustrated by Fig. 1, where 
the maps of the central part of the galaxy M82 are shown 
in both “high” and “low” states. 


3 DISCUSSION 
3.1 “Propeller” effect 

The remarkable behaviour of M82 X-2 showing dramatic 
switches in luminosity by a factor of 40 can be inter¬ 
preted as the onset of the so-called “propeller effect” 
(Illarionov & Sunyaev 1975). This effect is caused by a sub¬ 
stantial centrifugal barrier which have to be broken by the 
infalling matter during the accretion onto the rotating neu¬ 
tron star with strong magnetic field. At the magnetospheric 
radius where the magnetic pressure equals the ram pressure 
of the infalling material, the accreting matter from a disc or 
a wind is “frozen” into the stellar magnetic field lines and 
rotates rigidly with the angular velocity of the star. The 
matter will fall onto the neutron star only if the velocity of 
the magnetic field lines at the magnetospheric radius is lower 


^ The over-all bimodal flux distribution is confirmed by Fig. 3 
from Brightman et al. (2015). The only flux measurement there 
fallen between the two states is the Chandra observation 10545 
(MJD 55405), where the PSF of M82 X-2 is clearly blended with 
a nearby source, and hence is missing in our list of observations. 


than the local Keplerian velocity. Otherwise, the matter will 
be stopped at the radius of magnetosphere or even expelled 
from the system. Given the fact that magnetospheric radius 
depends only on the mass accretion rate and the strength 
of the magnetic field, the latter can be directly estimated if 
the propeller regime is observed in an accreting magnetized 
neutron star. 

The threshold value of accretion luminosity Luui{R) for 
the onset of the propeller can be estimated by equating the 
co-rotation radius (where a Keplerian orbit co-rotates with 
the neutron star) 


fCMP^Y^^ 
" ^ V 47r2 ) 


( 1 ) 


to the magnetospheric radius 

= fcM“^/'^//'^(2GM)“^/'^. (2) 


Here M is the neutron star mass, P its rotational period, /r = 
BR? 12 is the magnetic dipole moment related to the surface 
polar dipole magnetic field strength B and the neutron star 
radius R, and M is the mass accretion rate onto the neutron 
star. The dimensionless factor k relates the magnetospheric 
radius to the Alfven radius computed for spherical accretion; 
for disc accretion it is usually taken k = 0.5 (Ghosh & Lamb 
1978). At the limiting accretion rate M = Mum, Rc = Rm, 
and the accretion luminosity is (Gampana et al. 2002) 

Lnm(R) ^ ~ s"!, 

H 

( 3 ) 

where Mi .4 is the neutron star mass in units of I. 4 M 0 , Re is 
neutron star radius in units of 10® cm, B \2 is the magnetic 
field strength in units of 10^^ G. 

The decrease of the accretion rate below Mum will lead 
to the propeller regime of accretion. The accretion efficiency 
drops significantly and the luminosity in that regime cor¬ 
responds to the accretion onto the magnetosphere with the 
maximum value of (Corbet 1996) 

Liim(i?o) = ^ Lnm{R)^. (4) 

-ttc -^C 

Thus if the pulsar is close to the spin-equilibrium, when 
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Figure 2. (a) Light curve of M82 X-2 obtained by the Chandra observatory during 15 years of observations. Luminosities are corrected 
for the absorption and given in the energy range 0.5-10 keV; (b) Distribution of individual observations over luminosities (black line). 
Bimodal structure is clearly seen. Red and blue dashed lines show the averaged luminosities in the “high” and “low” states, respectively. 
Grey line represents rescaled distribution of luminosities of X-ray pulsar LMC X-4 from the Swift/HAT data; (c) Predicted dependence 
of the magnetised neutron star luminosity on the mass accretion rate for the following parameters; spin period P = 1.37 s, magnetic field 
strength B = 1.1 X 10^^ G, neutron star mass M = 1.4M0 and radius ij = 10 km, magnetospheric radius in units of the Alfvenic one 
k = 0.5. A small fraction (2.5 per cent) of the accreting material is assumed to leak through the magnetosphere onto the neutron star 
surface. 


Rc ~ Rm, small variations in the accretion rate will lead to 
large variations in the observed luminosity: 


LlimiR) _ 

LnmiRc) ^ 'v47r2i?3 




( 5 ) 


Exactly such behaviour, i.e. abrupt switches between 
two intensity states, is observed in M82 X-2 (see Fig. 2(b)). 
To our knowledge none of the known other variability mech¬ 
anisms would result in such a sharp bimodal distribution 
of flux from an accreting magnetized neutron star. For a 
comparison, we show in Fig. 2(b) by grey line the lumi¬ 
nosity distribution for the X-ray pulsar LMC X-4 demon¬ 
strating a well established super-orbital variability as seen 
by the SwiftfBAT all-sky monitor.® The variability patterns 
obviously differ much. This is also true for other sources 
with periodic (SMC X-1, Her X-1) and non-periodic (e.g., 
Cen X-3) variability of the flux. Unfortunately, the relatively 
short exposures of the Chandra observations do not permit 
us to detect transitions from one state to another. How¬ 
ever, a quite dense observational program of M82 allows us 
to put an upper limit on the duration of such transitions: 
two observations separated by 37.3 days (ObsID 1302 and 
1411_000) show flux change by a factor of ~ 50. 


3.2 Magnetic field measurement 

Theoretically, expressions (3)-(5) would allow us to measure 
the magnetic field of the neutron star together with its com- 

® See http: //swift .gsfc .nasa.gov/results/traiisients/LMCX-4/ . 
Daily averaged count rates were translated into the 
source luminosity using its broadband spectrum (see, e.g., 
Tsygankov & Lutovinov 2005). The obtained luminosity distri¬ 
bution was shifted by a factor of 20 to visually match the M82 
X-2 luminosity range. 


pactness. However, the source does not follow the theory so 
closely. Indeed, the observed ratio of luminosities in “high” 
and “low” states is 40, that is smaller than the ratio around 
210 predicted by equation (5) for a given pulse period, prob¬ 
ably because of a substantial leakage of accreting material 
through the magnetosphere. Such a behaviour is typical for 
accreting X-ray pulsars (Doroshenko et al. 2011, 2014). This 
is also supported by the fact that we do not detect signif¬ 
icant changes of the source spectra in the two states. The 
jump by a factor of 40 in the luminosity can be easily ex¬ 
plained if 2-3 per cent of accreting matter leaks through the 
magnetosphere (see Fig. 2c). 

However, the limiting luminosity in the high state is less 
affected by such uncertainties in the accretion model and can 
be used to estimate the magnetic field strength of the neu¬ 
tron star directly from equation (3). For this purpose we take 
P — 1.37 s and Lum = 2.0 x 10"^° erg s“® (where we apply a 
bolometrlc correction of factor 2 to the observed mean lumi¬ 
nosity in the high state). This gives immediately a magnetar- 
like magnetic field of B = (1.1T0.4) x 10®“^ G, which exceeds 
by a factor of 15 the largest measured so far magnetic field 
in an accreting X-ray pulsar (Yamamoto et al. 2013). The 
error on B comes from the uncertainty in a combination of 
factors 4 ^^Re in relation (3) as well as from a 50 per 

cent uncertainty in determination of the limiting luminosity 
(as a consequence of a finite width of the luminosity distri¬ 
bution in “high” state) and the bolometric correction. For 
example, if instead of using an averaged luminosity in high 
state as the limiting one we take the minimal observed in 
this state (1.2 x lO'^® erg s~®), the magnetic field strength 
would be B = 8.5 x 10^® G. 

The reliability of this method is illustrated by Fig. 3, 
where we show the correlation between the combination 
BiimB^^® and the magnetic field strength for one accreting 
millisecond pulsar SAX J1808.4-3658 and three X-ray pul- 
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Figure 3. The observed correlation between the magnetic field strength B and a combination of the limiting luminosity at the start 
of the propeller regime and the period, for four pulsars is shown by circles with error bars. The dashed line represents the 

theoretical dependence from equation (3) assuming standard parameters M = IAMq, i? = 10 km, k = 0.5. The star indicates the position 
of the ULX M82 X-2 for which we estimate = 5 X 10^° erg s"^/®, that corresponds to the B-field strength of 1.2 X 10^"^ G. 


sars GRO J1744-28, 4U 0115-1-63, and V 0332-1-53 where 
the action of the propeller was mentioned in the literature 
(see Appendix A). This sample consists only of sources with 
the confidently determined nature of the compact object as 
a neutron star with a sufficiently strong magnetic field, i.e. 
sources with the pronounced X-ray pulsations. We see that 
the data spread over four orders of magnitude in B are well 
described by a theoretical dependence given by equation (3) 
and shown by the dashed line in Fig. 3. 


4 CONCLUSION 

In spite of a number of suggestions for the presence of mag- 
netars in binary systems (Shakura 1975; Bozzo et al. 2008; 
Doroshenko et al. 2010; Reig et al. 2012; Klus et al. 2014), 
no such systems up to date have been unambiguously iden¬ 
tified. Here we show that the ULX M82 X-2 regularly enters 
the “propeller regime” of accretion which is seen as dramatic 
variations of the emitted luminosity and two well-defined 
peaks separated by a factor of 40 in the X-ray luminosity 
distribution. These observations imply the dipole compo¬ 
nent of the neutron star magnetic field of ~ 10^“* G, making 
the source the first confirmed accreting magnetar. 

Our discovery of such a strong magnetic field in a 
pulsar-ULX naturally explains the observed deficit of ULXs 
powered by accretion onto a neutron star. For instance, 
binary population synthesis models suggest (Shao & Li 
2015) neutron star-ULX to be not less numerous than 
black hole-ULX. Moreover, the observed binary parameters 
(Bachetti et al. 2014) of M82 X-2 (donor mass Me > 5.2 Mq 
and orbital period ~ 2.5 d) are shown to be typical. However, 
these studies do not account for an extreme narrowness of 
all physical parameters permitting for such a binary to oper¬ 
ate as ULX (Mushtukov et al. 2015). Namely, the magnetic 


field strength should be high enough to make the scattering 
cross-section small in order to support super-Eddington lu¬ 
minosity from the accretion column. Furthermore, for a high 
accretion rate needed for a neutron star to become ULX, 
the magnetic field should be high enough to make magne¬ 
tosphere larger than the spherization radius (where the ac¬ 
cretion disc scale-height becomes comparable to the distance 
from the neutron star), otherwise the accreting material will 
be blown away by radiation and the accretion may proceed 
only in a nearly spherically-symmetric fashion, which is lim¬ 
ited by the Eddington luminosity (Mushtukov et al. 2015). 
This naturally limits the magnetic field of pulsar-ULXs to 
be greater than ~ 3 x 10*^® G. 
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APPENDIX A: OBSERVATIONAL EVIDENCE 
OF THE ONSET OF THE PROPELLER REGIME 

For Fig. 3 we have collected all the cases of X-ray pulsars 
where the action of the propeller was mentioned in the lit¬ 
erature. We selected only sources with the confidently de¬ 
termined nature of the compact object as a neutron star 
with a sufficiently strong magnetic field, i.e. sources with 
the pronounced X-ray pulsations. 

The first discovered accreting millisecond pulsar SAX 
J1808.4-3658 has the shortest spin period of P = 2.5 ms 
and the lowest magnetic field strength among sources in 
the sample determined to be B = (0.8 ± 0.5) x 10*^“^^"* G 


(Ibragimov & Poutanen 2009). Substituting k — 0.5 (as we 
do in our work) the final magnetic field strength on the neu¬ 
tron star shown in Fig. 3 is P = (2.7 ± 1.7) x 10® G. Vari¬ 
ations in the flux by a factor of 100 have been observed 
by the Swift observatory during the 2005 outburst and in¬ 
terpreted as the transition to the propeller regime of ac¬ 
cretion (Gampana et al. 2008). The luminosity at the on¬ 
set of propeller in the 0.3-10 keV energy range was esti¬ 
mated to be L = (3 ± 1) X 10®®erg s“® assuming a source 
distance of 3.5 kpc (Galloway & Gumming 2006). Apply¬ 
ing a bolometric correction of factor 1.7 we get Lum = 
(5 ±2) X 10®® erg s"\ 

An intermediate X-ray pulsar GRO J1744-28 situated 
1° from the Galactic centre has the pulse period of 0.467 s 
and the magnetic field B ~ 5.3 x lO*^® G determined from 
the cyclotron absorption line at Pcyc ~ 4.7 keV (D’Ai et al. 
2015; Doroshenko et al. 2015). Because this feature origi¬ 
nates from the vicinity of the neutron star surface, its energy 
was corrected for the gravitational redshift to be compared 
to the magnetic field of SAX J1808.4-3658. The disappear¬ 
ing of strong pulsations from the source during the decay 
of 1996 outburst when the 2-60 keV flux was in the range 
F = (4±2) X 10~® erg cm~^ s“® was interpreted as the onset 
of the propeller regime (Gui 1997). Assuming a distance to 
the source of 8 kpc (Giles et al. 1996), the propeller thresh¬ 
old luminosity can be estimated as Lum = (3.0 ± 1.5) x 10®®" 
erg s“®. 

The propeller effect was observed also in two well- 
studied transient X-ray pulsars with Be companions - 
4U 0115-1-63 and V 0332-1-53, characterized by the spin pe¬ 
riods of 3.6 s and 4.35 s and the magnetic fields of 
B « 1.3 X 10®^ G (White et al. 1983) and 3.0 x 10®^ G 
(Makishima et al. 1990), respectively. The sudden decrease 
of the flux of 4U 0115-1-63 and its pulsed fraction below 
the luminosity Liim = (1 ± 0.5) x 10®® erg s“® can be 
interpreted as the onset of the propeller (Campana et al. 
2001). In the case of V0332+53, the minimum flux in 
the 1-15 keV range measured just before the turnoff was 
B = (3 ± 1) X 10“®® erg cm“^ s“® (Stella et al. 1986). Tak¬ 
ing the distance to the source of 7 kpc (Negueruela et al. 
1999) and applying the bolometric correction of factor 1.5 
we get Biim = (2.6 ± 0.9) x 10®® erg s“®. We use a 10 per 
cent uncertainty in the magnetic field strength for the three 
pulsars, because of the uncertainty in the redshift correction 
and in the measured energy of the cyclotron line. 

This paper has been typeset from a J+X/IAIliX hi® prepared by 
the author. 
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